The vibrational stability and electronic structure of B 8 o fullerene-like cage 
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We investigate the vibrational stability and the electronic structure of the proposed icosahedral 
fullerene-like cage structure of Bgo [Szwacki, Sadrzadeh, and Yakobson, Phys. Rev. Lett. 98, 166804 
(2007)] by an all electron density functional theory using polarized Gaussian basis functions con- 
taining 41 basis functions per atom. The vibrational analysis of Bgo indicates that the icosahedral 
structure is vibrationally unstable with 7 imaginary frequencies. The equilibrium structure has T; t 
symmetry and a smaller gap of 0.96 eV between the highest occupied and lowest unoccupied molec- 
ular orbital energy levels compared to the icosahedral structure. The static dipole polarizability of 
Bso cage is 149 A 3 and the first ionization energy is 6.4 eV. The Bgo cage has rather large electron 
affinity of 3 eV making it useful candidate as electron acceptor if it is synthesized. The infra-red 
and Raman spectra of the highly symmetric structure are characterized by a few absorption peaks. 
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Since the discovery of Ceo, several studies have re- 
ported possible existence of hollow inorganic cage-like 
structures. Recent work by Szwacki etal[l| added boron 
to the list of elements that can form fullerene-like hollow 
cag e structure. While several studies @, 0, 0, 0, 0, H, @, 
lfj have reported stable clusters, rings, and nanotubes of 
boron, a hollow fullerene-like cage cluster containing only 
boron has not yet been found. Using density functional 
theory Szwacki etal[l[, showed that a boron cluster con- 
taining 80 atoms can form a stable hollow cage. The basic 
structure of this cluster is similar to that of Ceo with 12 
pentagonal and 20 hexagonal rings. The Bgo cluster has 
an additional boron atom at the center of each hexagon 
(Cf. [T]). It is known that the boron analogue of Ceo is not 
a stable structure. However, these additional 20 boron 
atoms at the centers of the hexagonal rings stabilize the 
Bgo to form a stable Bgo fullcrcnc with a binding energy 
of 5.77 eV per atom. The structure of Bgo cluster has 
been suggesteded by Szwacki etal. to maintain the same 
icosahedral point group symmetry as in Cgo fullerene. 
They also predicted that the icosahedral Bgo, as in the 
case of Ceo, has the lowest unoccupied molecular orbital 
(LUMO) of ti u symmetry. However they find that the 
highest molecular orbital (HOMO) however belongs to 
the t2u irreducible representation[l|. 

The purpose of the present article is to further inves- 
tigate the electronic structure of this novel boron clus- 
ter, focusing particularly on its response to applied elec- 
tric field. To this end, we calculate the dipole polariz- 
ability of the Bgo cluster using the linear combination 
of atomic orbital approach within the density functional 
formalism. In addition to the dipole polarizability we 
also compute the first ionization energies and the elec- 
tron affinity. We also provide the infra-red and Raman 
spectra for possible detection. We use large polarized 




FIG. 1: (Color online) Optimized geometry of Bgo cage. 



Gaussian basis sets to express the Kohn-Sham molecular 
orbitals[Hj]. The exchange-correlation effects are treated 
within the generalized gradient approximation using the 
Perdew-Burke-Ernzerhof[12j parametrization. We have 
first optimized the structure of B 80 using the icosahedral 



point group by the NRLMOL code[13lll4|. The resultant 
electronic structure shows some discrepancy with respect 
to the previous study by Szwacki [lj]. The HOMO level is 
5 fold degenerate and belongs to the h u irreducible repre- 
sentation instead of t2« as reported earlier. Besides, this 
small discrepancy, the present calculation also revealed 
that a more stable Bgo cluster has symmetry lower than 
the icosahedral symmetry as reported by Swazcki et oi. . 
To verify if this lower symmetry structure is due to sym- 
metry breaking distortion of icosahedral cage, full vibra- 
tional frequency calculations were performed within the 
harmonic approximation. These calculations show that 
the icosahedral Bgo cluster is vibrationally unstable with 
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FIG. 2: (Color online) Unstable vibrational modes of icosa- 
hedral Bso cage. 



two sets of imaginary frequencies. One of the two sets of 
modes is four fold degenerate and the the other one is a 
three fold degenerate mode. The vibrational motions of 
atoms in these two modes are depicted in Fig. [2] These 
modes are the out of plane motion of the boron atoms 
capping hexagons. 

Full symmetry unconstrained relaxation of Bgo cluster 
show that the cluster has Th symmetry. The Th clus- 
ter is vibrationally stable. The Th Bgo cluster has five 
inequivalent atoms whose positions in Bohr are (0.0000, 
7.9117, 1.6024), (4.3303, 4.3303, 4.3303), (2.6304, 6.8420, 
3.2061), (5.1888, -5.8587, -1.6209), and (0.0000, 2.5093, 
6.6176). The position of all atoms in Bgo cluster can be 
obtained using the position of the nonequivalent atoms 
and using the symmetry operations of Th- The 4 three 
fold symmetry axes are along the 111 directions. 

The energy of T^ structure is lower by 0.05 eV com- 
pared to the icosahedral structure. Its electronic struc- 
ture is 7& u 12& g 22e g 22t 9 27t u 12e„. The h„ symmetry 
of the HOMO of the icosahedral structure is split into 
a three-fold degenerate t u and a two fold degenerate e u 
level in the Th structure. The splitting also reduces the 
HOMO-LUMO gap in the T h structure. The HOMO is 
two fold degenerate and belongs to e u irreducible repre- 
sentation. The LUMO is three fold degenerate is of t u 
type. The HOMO and LUMO orbital densities are shown 
in FigHl 

The static dipole polarizability is an important phys- 
ical property that characterizes the system's response 
to an applied static electric field to the first order. A 
number of methods have been developed to compute 
the dipole and higher polarizabilities. In this work, we 
use the finite-field method. In the finite-field method 
the total energies and/or dipole moments are computed 
self-consistently for various values of the applied electric 
field [H [H, [H 0, H, Ei The polarizability tensor is 
then built from well converged total energies or dipole 
moments using the finite-difference approximation. An 
accurate estimate of the polarizability using the linear 
combination of atomic orbitals requires the use of a large 
basis set supplemented with diffuse functions. The Gaus- 



sian basis set used in this work consists of 5 s— , 4 p— , 
and 3 d— type Gaussians each contracted from 12 prim- 
itive functions. This basis is augmented with 6 d— type 
functions. Thus, in total 3280 basis functions are used in 
polarizability calculation. More details about construc- 
tion and theperformance of basis set can be found in 
Ref. [ll|, [U [22]. Due to the quasi spherical symmetry of 
Bgo cluster, the off-diagonal elements of the polarizabil- 
ity tensor are zero. The mean polarizability obtained by 
finite-field method is 149 A 3 . The mean polarizability of 
Ceo fullerene determined using the same set of approxi- 
mations is 82 A 3 . The larger polarizability of Bgo cluster 
is principally due to its larger volume. Using classical 
electrostatics, the polarizability of a spherical shell of ra- 
dius R can be shown to be R 3 . Unlike the Ceo fullerene 
where in all atoms are at same distance from its center 
of mass, the atoms in Bgo cluster are at slightly different 
distances (3.74-4.27A) from the center of Bgo cluster. Us- 
ing the radius of the outermost atoms, the volume of Bgo 
cluster can be estimated to be roughly 1.74 times that 
of the Ceo fullerene. Using 1.74 as a scaling factor, the 
polarizability of Cgo fullerene can be used to estimate po- 
larizability of Bgo cluster. This rough estimate is 143 A 3 , 
in good agreement with the polarizability obtained by 
finite-field method. The finite-field (screened) polariz- 
ability does not change due to symmetry lowering of the 
Bgo cluster. The unscreend polarizability however show 
significant change upon symmetry lowering. It increases 
from 63lA 3 to 897 A 3 due to the decrease in the HOMO- 
LUMO gap and due to changes in low-lying dipole al- 
lowed transitions. These transitions must be strongly 
screened to give identical values of polarizability for Ih 
and Th clusters. 

The calculation of the vibrational frequencies estab- 
lishes the stability of the T^ structure. The vibrational 
density of states of the Bgo cluster are shown in fig. [3] 
The bottom panel shows the density of states and the 
upper two panels show respectively the infrared and Ra- 
man activity of the Bgo cluster. The cluster shows a very 
strong infrared peak at 991 cm . The other significant 
peaks occur at 446, 759, 771, 912, and 1012 cm- 1 . The 
Raman spectrum shows a few low frequency modes with 
strong peaks at 112, 174, 176, and 312 cm" 1 -1 . An- 
other prominent peak is seen at 963 cm . 

The vibrational frequencies determined in analysis of 
the vibrational stability of Bgo cluster can be also be used 
to compute the vibrational contribution to the dipole po- 
larizability. In general, the vibrational polarizability is 
the second largest contribution to the total polarizability 
tensor. For the case of ionic and hydrogen bonded sys- 
tems, the vibrational contribution may be comparable to 
or even larger than the electronic contribution. The vi- 
brational contribution to the polarizability tensor within 
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FIG. 3: (Color online) The vibrational density of states, In- 
frared absorption intensity and the Raman scattering inten- 
sities of the Bgo cluster with structure. 

the double harmonic approximation [22} is given as 

Vib _ \ "* ry -2 ryT 

Here, u>~ 2 is the frequency of the /xth vibrational mode, 
Zi 4L is the effective charge tensor (See. Ref. [22] for more 
details) . The vibrational contribution to polarizability is 
5.5 A 3 which is much smaller the electronic contribution 
but is larger than that observed in carbon fullerenesf^. 

To understand the reactivity of the Bso cluster, we 
have calculated its ionization potential and the electron 
affinity. The first ionization energy is the minimum en- 
ergy required to remove an electron. It can be deter- 
mined from the self-consistent total energy difference of 
the Bgo cluster and its singly charged cation[23]. The 
calculated ionization energy is 6.5 eV and is smaller than 
that of Cgo (7.6eV). Similarly, the electron affinity can 
be computed from self-consistent total energy difference 
of neutral cluster and its anion. The calculated elec- 
tron affinity is 3 eV. Our calculations also indicated that 
the cluster reorganization upon addition of an electron is 




FIG. 4: (Color online) The orbital densities of the HOMO 
(left) and LUMO (right) of the B g0 cluster. 



also small in this highly symmetric cluster. The chemical 
hardness is an indicator of the reactivity of the molecule. 
It can be approximated as half of the difference between 
the ionization potential and electron affinity [24]]. This 
definition immediately points out that the Bso cluster is 
more reactive than the Ceo fullerene. The rather large 
electron affinity makes B§o cluster, if synthesized, an in- 
teresting candidate as an electron receptor. The larger 
electron affinity also suggests possibility of coating these 
cluster with alkali or transition atoms, which then can be 
tested as hydrogen storage materials. Other possibilities 
are using it as a building block in ionic cluster assembled 
materials or as an electron receptor in a photovoltaic de- 
vice. 

To summarize, the vibrational stability of recently re- 
ported Bso cluster is examined by computing the har- 
monic vibrational frequencies. The Bso cluster is found 
to be unstable in icosahedral symmetry but stable in the 
reduced Th symmetry. The electronic structure of Bgo 
cluster changes upon symmetry lowering. The symme- 
try lowering decreases the HOMO-LUMO gap from 1.1 
eV (in I/j structure) to 0.97 eV. Its static dipole polariz- 
ability (149A 3 ) however does not alter appreciably. Its 
ionization potential is 6.5 eV. The structure has large 
electron affinity of 3 eV making it a candidate as an elec- 
tron receptor. Using the chemical hardness as an indica- 
tor of reactivity the Bso is found to be larger than the 
Ceo fullerene. The infra-red and Raman spectra are pro- 
vided. 

This work is supported in part by the National Sci- 
ence Foundation through CREST grant, by the Univer- 
sity of Texas at El Paso (UTEP startup funds, University 
research institute grant) and partly the Office of Naval 
Research, directly and through the Naval Research Lab- 
oratory Authors acknowledge the computer time at the 
UTEP Cray acquired using ONR 05PR07548-00 grant. 



[1] N. G. Szwacki, A. Sadrzadeh, and B. I. Yakobson, Phys- 
ical Review Letters 98, 166804 (2007), pT: J. 

[2] V. Dadashev, A. Gindulyte, W. N. Lipscomb, L. Massa, 
and R. Squire, Structures and Mechanisms: from Ashes 
to Enzymes 827, 79 (2002), pT: S. 

[3] A. Derecskeikovacs, B. I. Dunlap, W. N. Lipscomb, 
A. Lowrey, D. S. Marynick, and L. Massa, Inorganic 
chemistry 33, 5617 (1994), pT: J. 

[4] A. Gindulyte, N. Krishnamachari, W. N. Lipscomb, and 
L. Massa, Inorganic chemistry 37, 6546 (1998), pT: J. 

[5] A. Gindulyte, W. N. Lipscomb, and L. Massa, Inorganic 
chemistry 37, 6544 (1998), pT: J. 

[6] H. Tang and S. Ismail-Beigi, Physical Review Letters 99, 
115501 (2007), pT: J. 

[7] S. Chacko, D. G. Kanhere, and I. Boustani, Phys. Rev. 
B 68, 035414 (2003). 

[8] K. C. Lau, R. Pandey, R. Pati, and S. P. Kama, Ap- 
plied Physics Letters 88, 212111 (pages 3) (2006), URL 



4 



http : //link . aip . org/link/?APL/88/212111/l 

[9] I. Boustani, A. Quandt, E. Hernandez, and A. Rubio, [18] 
The Journal of Chemical Physics 110, 3176 (199 9), URL 
|http : //link . aip . org/link/? JCP/1 10/3176711 [19] 

[10] I. Boustani and A. Quandt, EPL (Eu- 

rophysics Letters) 39, 527 (1997 ), URL [20] 
|http : //stacks . iop . o rg/0295-5075/39/527] 

[11] D. Porezag and M. R. Pederson, Physical Review a 60, [21] 
2840 (1999), pT: J. 

[12] J. P. Perdew, K. Burke, and M. Ernzerhof, Physical Re- [22] 
view Letters 77, 3865 (1996), pT: J. 

[13] M. R. Pederson and K. A. Jackson, PHYSICAL RE- 
VIEW B. CONDENSED MATTER 43, 7312 (1991), 
pUBM: Print; ppublish. [23] 

[14] M. R. Pederson and K. A. Jackson, Physical Review B 

41, 7453 (1990), pT: J. [24] 

[15] S. A. Blunder!, C. Guet, and R. R. Zope, Phys. Rev. Lett. 
84, 4826 (2000). 

[16] K. Jackson, M. Pederson, C.-Z. Wang, and K.-M. Ho, 
Phys. Rev. A 59, 3685 (1999). 

[17] H. D. Cohen and C. C. J. Roothaan, The Journal of 



Chemical Physics 43, S34 (1965). 

R. J. Bartlett and G. D. Purvis, Phys. Rev. A 20, 1313 
(1979). 

H. A. Kurtz, J. J. P. Stewart, and K. M. Dieter, J. Corn- 
put. Chem. 11 (1982). 

J. Guan, M. E. Casida, A. M. Koster, and D. R. Salahub, 
Phys. Rev. B 52, 2184 (1995). 

M. R. Pederson and T. Baruah, Lecture Series in Com- 
puter and Computational Sciences 3, 156 (2005). 
M. Pederson, T. Baruah, P. Allen, and C. Schmidt, 
Journal of Chemical Theory and Computa- 
tion 1, 590 (2005), ISSN 1549-9618, URL 

http : //pubs3 . acs . org/ acs/journals/doilookup?in_doi=10 . 1021/ 
R. O. Jones and O. Gunnarsson, Reviews of Modern 
Physics 61, 689 (1989), pT: J. 

R. G. Parr and R. G. Pearson, Jour- 
nal of the American Chemical Society 

105, 7512 (1983), ISSN 0002-7863, URL 

http : //pubs3 . acs . org/ acs/journals/doilookup?in_doi=10 . 1021/ 



